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Abstract 

We calculate low scale gravity effects on the cross section for neutrino-nucleon scattering 
at center of mass energies up to the Greisen-Zatsepin-Kuzmin (GZK) scale, in the eikonal 
approximation. We compare the cases of an infinitely thin brane embedded in n — 5 compactified 
extra-dimensions, and of a brane with a physical tension Ms = 1 TeV and Ms = 10 TeV. The 
extra dimensional Planck scale Mjj is set at 10 3 GeV and 2 x 10 3 GeV. We also compare our 
calculations with neutral current standard model calculations in the same energy range, and 
compare the thin brane eikonal cross section to its saddle point approximation. New physics 
effects enhance the cross section by orders of magnitude on average. They are quite sensitive to 
Ms and Md choices, though much less sensitive to n. 

1 Introduction 

Since neutrinos interact only weakly with matter, neutrino observatories provide a powerful tool 
for exploring the deepest reaches of stars and galaxies. In recent years there has been an upswell in 
the number of experiments, either running or about to run, that are aimed at detecting ultra-high 
energy (UHE) neutrinos, neutrinos with energies from the multi-TeV range to beyond the EeV 
range. Their existence is predicted by a number of theoretical models. Although the production 
mechanism can be different - either as a byproduct of creation and decay of pions in the interaction 
between high energy primary cosmic rays and the cosmic microwave background, as predicted by 
GZK [lj, or as products of the same astrophysical sources that generate the cosmic rays observed 
in the highest energy air showers [2] -there is wide expectation in the community that neutrinos 
of galactic or cosmic origin will be detected in the UHE energy range in the near future. Many 
experimental limits on astrophysical neutrino fluxes have already been established. Frejus, Baikal, 
AMANDA and MACRO [3] have reported limits on neutrinos from astrophysical sources in the 
10 3 - 10 6 GeV range, while Fly's Eye, AGASA, AMANDA and RICE [4j have given limits in the 
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range 10 6 - 10 9 GeV, and AGASA, Fly's Eye, RICE, Forte, ANITA and GLUE [5] have reported 
limits in the range above 10 9 GeV. In the near future, IceCube [6], ANITA and AUGER [7] can 
be expected to release stronger results based on more recent data. In the longer term, proposed 
expansions of IceCube, such as AURA [8j, or a salt-based radio telescope like SALSA [9] could afford 
substantially enhanced sensitivity. Further in the future, an orbiting telescope like the proposed 
EUSO project [10] opens the possibility to achieve a huge effective volume. 

The design of neutrino telescopes depends critically on the estimates of neutrino cross sections 
in the UHE regime, most of which is far beyond currently available data. Moreover, the systematic 
study of the scattering of such high energy neutrinos with baryonic matter might prove to be an 
essential instrument to test new physics effects, such as those given by models of low scale gravity 
(LSG). Models of extra dimensions \12\ [T3] have lowered the characteristic quantum-gravity 
scale to energies comparable with the electro-weak scale, enhancing the expected neutrino cross 
sections at GZK-energies and above, even for conservative parameter choices. 

Efforts at realistic calculations of the UHE neutrino-nucleon cross section have been performed in 
the framework of QCD and the standard model |14t ll5j. and many speculative calculations based on 
new physics effects have been presented over the past decade or more. The only detailed calculations 
of neutrino event rates within the framework of the LSG models that include the neutrino-proton 
eikonal cross section with LSG graviton exchange, rely upon the saddle point approximation |16j . 
In this article we present in detail the calculation for the full eikonal approximation of the LSG 
neutrino-proton cross section, including a study of the effects of a finite tension on the brane [171118]. 
We comment on the inapplicability of the saddle point approximation at neutrino energies below 
10 11 GeV. 

In Sec. II we summarize the main procedure of the eikonal approximation (EA) in quantum field 
theory and apply it to the calculation of the desired differential cross sections at the parton level 
in the Arkani-Hamed, Dimopoulous and Dvali (ADD) model of extra dimensions [11]. We consider 
the cases of an infinitely thin standard model brane and the one where it has finite thickness. In 
Sec. Ill we present the full EA calculation and compare the results with those in the saddle point 
approximation, and with the standard model cross section. Sec. IV is dedicated to the summary of 
results, comments on the sensitivity of results to choice of gravity scale and brane tension and on 
the origin of the discrepancy between the saddle point and full EA calculation. We also comment on 
the impact of our cross section results on neutrino telescope event rate expectations. The Appendix 
presents details of the representation of the eikonal amplitude in terms of Meijer G-functions. 
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2 Theory 

The eikonal approximation is a technique widely used in quantum mechanics and wave physics to 
derive an expression for elastic scattering in the limit of large center of mass (CM) energies (see for 
example Reference |19j). In quantum field theory the amplitude can be obtained by summation of 
all the ladder and cross-ladder diagrams for boson exchange at all order [20\ [2~T] , as exemplified in 
Fig. 1. 




Figure 1: Sum of the ladder and cross-ladder diagrams in the eikonal approximation 

The approximation is valid for center of momentum energies that are high with respect to Q 2 , 
the absolute magnitude of the momentum transfer q 2 , i.e. \t/s\ < 1 in the usual Mandelstam 
variables for the 2^2 scattering. The high energy scattering angle is small and the four-momenta 
of the incoming particles are approximately equal to those of the corresponding outgoing particles, 
pi ~ p'i and p2 ~ P2- For any given kind of boson exchanged (scalar, vector, tensor), the vertex 
and external legs factor is independent from the spin of the incoming particles, as one can neglect 
their mass and the recoil of the matter field. 

With these assumptions, s ~ 2(pi ■ P2) and the tree-level amplitude (Fig. 2a), or Born term for 
the exchange of one intermediate boson is given by 

I 

A B orn(t) = g 2 S r - 5-, (1) 

t — m z 

where g is the dimension-dependent interaction coupling constant, m is the mass of the exchanged 
particle and r = 0, 1, 2 is its spin. 

At one-loop level (Fig. 2b), the EA simplifies the matter propagators: 



(p ± q') 2 -M 2 + ie ±2p -q' + ie' w 
where M is the mass of the incoming particles. The one-loop eikonal amplitude, obtained when 
the above propagators are put on shell by the Cutkowsky rule, is imaginary. The integral over the 
transferred four-momentum reduces to an integral over its components in the plane perpendicular 
to the momenta of the incoming particles, q±. One gets for the eikonal at one loop: 
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(a) 





Figure 2: 2a: Born term diagram. 2b: One loop diagram with matter propagators put on shell by 
the Cutkowsky rule. 2c: Example of an inelastic diagram. 



= i/ 7^Jp" A ^ ow («'x)^ 4 s rn [(g_L - </l) 2 ] 

= -2sz | d 2 b ± e^ (-^x 2 ) (3) 
where the eikonal phase 

= ^ / ^e^W,?) (4) 

is the Fourier transform of the Born term. 

The procedure can be repeated with higher order diagrams [2D], such that one can sum all the 
terms in \ to obtain : 



A eik = Abotu + A loop + ... + A n _ loops + ... = -2si / d 2 b ± e iq ^ ( ^(ix) n /n! 




(5) 

It is worth noting at this point of our discussion that the validity of the EA in describing the 
asymptotic behavior of elastic scattering has been long debated. References [20] and [22] argue 
that the eikonal is not a correct description of the high energy limit in the cases of scalar particle 
exchange, abelian vector exchange and non-abelian vector exchange. In the cases of vector particle 
exchange, when we put r = 1 in Eq. ([T]), the EA describes accurately the asymptotic behavior of 
the sum of exchange diagrams, but fails to consider the inelastic diagrams related by unitarity to 
particle production, of which Fig. 2c is an example. As has been shown by [20], such diagrams 
become dominant at eighth order in g. On the other hand, the s 2 dependence of the graviton- 
exchange Born term in the quantum gravity EA renders the sum of exchange diagrams dominant 
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with respect to the inelastic unitarity diagrams ([22] and references therein). Consistent with the 
state of present knowledge, we will consider the eikonal a good approximation of the asymptotic 
behavior for graviton exchange. 

In models of LSG with large extra dimensions |11[ I13j. graviton exchange (r = 2) is ultraviolet 
divergent already at tree- level |23[ [23] . This is due to the presence of an infinite tower of massive 
Kaluza-Klein modes over which the amplitude has to be summed. For a number n of extra- 
dimensions we have, approximating the sum by an integral, 



A 



Born 



M™+ 2 



Sr, 



00 m n - x dm 



(6) 



t-m 2 M™ +2 Jq t-m 2 ' 
where Mp is the D-dimensional Plank mass (where D = 4 + n) and S n is the surface of a n- 
dimensional unit sphere. As the notation indicates, this sum is adopted as the "Born input" to the 
eikonal calculation |25j . 

In the literature, there are two common solutions to the problem of tree-level divergencies. One 
can introduce a cut-off of the order of the low gravity scale, as this is the limit below which the 
effective theory is deemed to be valid [231 12H [26]. Or, one can consider the physically reasonable 
case where the brane has a finite tension Ms, corresponding to an extension 1/Mg of the standard 
model fields along the extra dimensions [17[ I18j. In this case, if the extra-dimensional part y of the 
standard model wave function has a gaussian cut-off of the kind 



Eq. 



is modified to 



2vr/ 



(7) 



The Born term in Eq. 
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m n 1 dm 



io t — m l 

can be calculated by dimensional regularization 



(8) 



A 



Born 



^r(i-f) 



Substitution of the latter into Eq. (jlj) yields 



(9. 



2 x a _i 2 

^2 -V. 



0) 



X 



7T2 - x r(l - n/2)s 



4M™+ 2 



n-l 



Jo(qb) 



(4*)* _1 *r(3; 



2M™+ 2 



1 



b, 



(10) 



b n \b 

where we've followed the commonly simplified notation q = q±, b = b±. The critical impact 
parameter b c separates the region of space where the Born term is dominant (b 3> b c ) from the one 
where higher order terms in \ become relevant. Expression (|10p is clearly finite for all b > 0. Thus, 
it can be seen from Eqs. (f6|)- (|10j) that when the brane is infinitely thin (Ms — * oo) the ultraviolet 
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divergence gets absorbed by the eikonal phase calculation. Each regularization procedure leads to 

the same finite result [26]. 

The eikonal amplitude is then given by 



A eik = 4nsb 2 c F n (b c q), 



with 



(11) 
(12) 



F n (r,)=-i / d^J (^)(e^ -1), 
Jo 

and Jo the Bessel function of order zero. 

As mentioned in passing in [261 127] . but not pursued, the F n can be expressed in terms of Meijer 
G-functions. We implement our calculation explicitly in these terms, providing details in the 
Appendix. The result is 

F n (rj) = 2-%- 1 ir 1 2n- 1 (R n (r ] ) + il n (rj)), (13) 
where the real and imaginary part are given by 

„2n 

1 



T („\ — r<n+l,0 



and 



7] 



2n 



(n+1) I 22n+2 n 2n 
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(14) 



(15) 



where G^j are the Meijer G-functions |28| . 

In the limit of large momentum transfer with respect to the inverse impact parameter (q S> 
l/b c ), the amplitude of Eq. (jlip can be approximated analytically by the steepest descent method 
(PUGEl EZ]) and it reads 



i<t>n 



n+2 



qM D 



where 



3n 

(4vr) 2n+2 

y/n+1 



r. = + i 



1 

n + 1 



7T 



(>n = ^ + ( n + !) 



qb c 



n 



(16) 



(17) 



(18) 



and b c is given by Eq. (fTUj) . The location of the saddle point, b s =b c (-£-) n+1 , should satisfy b s <C b c 
in order for the saddle point to give the dominant contribution to the amplitude. 



On the other hand, in the case of a finite tension Ms in the brane, substitution of Eq. 
Eq. @ yields [H]: 

'b c M s ^ n 



sMa ^ /"n\ 7T2 



-1 



into 



(19) 
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where U(a, c; x) is a confluent hyper geometric function of the second kind [28J. 

By substitution of Eq. (fl~9j) into Eq. ([5]) the eikonal amplitude can be calculated easily: 



A e ik = — 47ris / dbbJo(qb) 
Jo 



be Ms ' 



(20) 



where in the last expression we have used the functional identity: 



2vr / dr}r}J (£rj)f(ri), 



(21) 



with £ = |£|, rj = \ ff\. We note here that, to our knowledge, there is no saddle point approximation 
for the extended brane case in the literature. Nor we were able to find for Eq. (|20p a simple, useful 
closed- form like (fl"6)) - (fl"8"|) . as the solution requires the numerical treatment of a transcendental 
equation. 

The differential cross section at the parton level can now be calculated with the usual substi- 
tution —t = —q 2 = Q 2 = xys, where x is the fraction of momentum of the proton carried by the 
parton and y is the inelasticity. We get for the infinitely thin brane: 

d 2 a 



dxdy 



^xfi(x,xys)snbc(s)\F n (b c (s)^/xljs)\' 



(22) 



with F n (rj) given by Eq. (fT2 
For the brane 1 /Ms-thick: 



d 2 a 
dxdy 



xfi(x, xys)s7r\G n (^/xys)\ 2 , 



with 



G n (0 



dbbJ {b() 



(23) 



(24) 



The fi(x,Q 2 ) are the usual parton distribution functions (PDF) summed over all the quark and 
antiquark flavors and the gluons. Note here that the parameter b c (s, n) of Eq. (|10p has to be 
calculated with respect to the center of mass energy squared of the neutrino-parton system s = xs. 



3 Results and Discussion 

The total cross section is obtained by integration of Eqs. (I22p and (|23f) . The y-integral has been 
performed with Mathematica using the CTEQ5 PDFs [29J, in the interval [Qq/(xs), l/(xsRg)], 
where Q 2 Q = O.OIM^, M w is the mass of the VF-boson and R% oc s V("+i)/ A fJ (n+2)/(n+1) is the 
Schwarzschild radius for black-hole production in the n-dimensional theory |26j . We adopt \/Q 2 
as the scale parameter in the PDF's. The upper limit of integration is introduced since the eikonal 
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cross section is not applicable for Q 2 > The lower limit, as in standard model calculations, 

allows us to be safely inside the perturbative QDC regime. The x-integral has been performed by 
dividing the interval [Qg/s, 1] into 38 logarithmic subintervals and performing Simpson's rule. This 
procedure gives an accurate evaluation of the cross section, which is reasonably smooth in x. 

We consider the case of n = 5 extra-dimensions, which is sufficient to illustrate all our points, 
since by changing n the calculations differ by amounts small compared to changes due to Md, M$- 
In Fig. 3 we show the comparison between our calculation of the proton-neutrino cross section in 
the full EA with the saddle point approximation of the same cross-section given by Eqs. (I16p . (|17j) 
and (|18|1 ^. The brane is here assumed to be infinitely thin, with the (4+5)-dimensional Planck 
scale set at Mb = 1 TeV. In the range E u > 10 7 GeV, where the EA is expected to be valid, 
the full eikonal cross section is significantly larger than its evaluation at the saddle point. The 
difference ranges from a factor 2.6 at the low end to about 1.5 at the high end. As discussed in 
the next section, the conditions for the dominance of the saddle point evaluation are not satisfied 
until extremely high energies, explaining the slow convergence of the saddle point value to the full 
eikonal result. We point out here that the calculation of the saddle point approximation at n = 6 
differs from the one at n = 5 by 10 — 20% over the entire energy range, in accordance with the 
findings of Ref. [16]. As it will be even clearer in what follows, this is a small amount in comparison 
to how the other parameters affect the results. 
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Figure 3: Solid light line: Full EA corrected by the Born term 1 , = 10 3 GeV, Ms = oo, 
n = 5. Dashed light line: Saddle point approximation corrected by the Born term. Solid dark line: 
Standard model neutral current. 

1 For an appropriate comparison with the calculations of Ref. [16], when evaluating the total cross section wc 
replaced the eikonal amplitude in Eq. ((22} by the Born term for x < Mp/s. For a detailed treatment of the Born 
term, see the first paper of Ref. |16j . 
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As also shown in Fig. 3, even with the saddle point approximation to the LSG eikonal, the 
neutral current cross section in this model exceeds the standard model values [15] at CM energies 
above the gravity scale, corresponding to roughly 10 6 GeV for Mp = 1 TeV [16J. The dominance 
is even more striking when the full eikonal is used: at neutrino energies greater than 10 10 GeV new 
physics effects raise the cross section by almost three orders of magnitude. 

The presence of a thick brane with a tension Ms = 1 TeV sharply reduces the cross section, 
due to the effective gaussian cut-off in the Born amplitude of Eq. ([8]) . The lower tension suppresses 
the eikonal cross section by an order of magnitude at the low energy end and a factor of two at the 
high end, as shown in Fig. 4. As the tension rises, and the brane gets thinner, the cross section 
should approach the infinitely thin limit. We already obtain an almost perfect overlapping when 
the tension is Ms = 10 TeV. 




10 7 10 s 10 9 10 10 10 11 

E, (GeV) 



Figure 4: Eikonal cross sections at LSG Mjy = 10 3 GeV for branes with different tension. (Not 
corrected by the Born term.) Solid light: Ms = 00, n = 5. Dashed light: Ms = 10 3 GeV, n = 5. 
Solid dark: Standard model NC. 

Fig. 5 shows the reduction to the cross sections due to a higher Planck scale M73 = 2 TeV. 
We here display the cases of the infinitely thin brane and of the Ms = 1 TeV-tension brane. We 
see that the increase in scale to 2 TeV suppresses the cross section an order of magnitude over the 
entire energy range for the thin brane case (solid light line). The relative suppression is even bigger 
in the case of a thick brane (dashed light line), where the cross section is lowered to values beneath 
the standard model calculation at energies less than 10 9 GeV. 
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Figure 5: Eikonal cross sections at LSG Mjj = 2 x 10 3 GeV for branes with different tension. (Not 
corrected by the Born term.) Solid light: Ms = do, n = 5. Dashed light: Ms = 10 3 GeV, n = 5. 
Solid dark: Standard model NC. 

4 Summary, Further Discussion and Conclusions 

We have numerically calculated the parton level cross-sections for UHE neutrino-nucleon gravi- 
tational scattering in the eikonal approximation, under the assumption that the standard model 
particles are confined in a 4-dimensional brane embedded in a higher number of compactified extra- 
dimensions (ADD Model). The analytical calculations are presented in the most general case of n 
extra-dimensions, but the numerical calculations were performed for n = 5 to illustrate the mag- 
nitude of the effects for a value of the extra dimensions that is only weakly constrained by current 
laboratory and astrophysical data (n > 5). Lower dimensions are required by astrophysical data to 
have corresponding gravity scales too high to allow significant effects on the UHE cross sections, 
whereas n > 5 are constrained only by collider data, which put lower limits on the gravity scale less 
than a TeV [301 13 lj . The cross sections are only mildly dependent on n, at any rate [16]. Previous 
applications of the ADD-type models to the UHE neutrino cross section and expectations for event 
rates relied on the saddle point estimation of the eikonal amplitude, and only in the thin brane 
limit. Our application completes the picture of the impact of low scale gravity on the predictions of 
high energy neutrino-nucleon cross sections. We reach the same qualitative conclusions as previous 
studies: the effects are sensitive to the scale, or radius, of the extra dimensions but much less 
sensitive to the actual number of dimensions. We also find that the brane tension, or thickness, 
greatly affects the size of LSG influence on UHE neutrino scattering, weakening sharply as brane 
tension decreases, or equivalently as brane thickness increases. This is in agreement with similar 
predictions for hadron collider applications |18j . 
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Quantitatively speaking, we have found that the presence of new physics at LSG M D = 10 3 
GeV enhances the "neutral current" cross section compared to the standard model cross section 
by three order of magnitude at CM energies of 10 10 GeV and above, and that the full eikonal 
calculation of the cross section gives values consistently higher than the corresponding saddle point 
approximation, when the brane is considered to be infinitely thin. The enhancement above the 
standard model expectation is not nearly as marked in the case of a brane physically extended 
10~ 3 GeV -1 along the extra-dimensions, or in the case where the LSG Planck Mass is increased 
by a factor of two. However, we found that calculations of the cross sections performed with brane 
tensions of 10 TeV and above are indistinguishable, indicating that a ratio of tension to gravity 
scale of 10:1 is effectively "infinite" for these purposes. 

The preceding statement is quite n-independent as can be seen from the r.h.s. of Eq. (|19p : since 
in the region of validity of the EA (s > M^) the impact parameter b is of the order of b c , one can 
see that when Ms 3> 6^7 1 the hypergeometric function U tends rapidly to (Msb c /2)~ n . The eikonal 
phase is therefore of order one, as expected, independently from the number of extra-dimensions. 
Correspondingly, the approach to the thin brane limit is only weakly dependent on n. 

To understand the difference between the UHE cross section evaluated in the saddle point 
approximation versus the full eikonal calculation, we recall that the saddle point dominance of the 
integral in question requires qb c 3> 1. In our eikonal evaluation, as in those in Ref. [IB], we required 
q < 1/Rs- Let us check the consistency of these two requirements. Defining q ma x = we find 

that 0.77 < q max x b c < 5.1 as E v runs from 10 5 GeV to 10 12 GeV for M D = 10 3 GeV and n = 5. 
The ranges are 0.71 < q ma x x b c < 8.0 and 0.81 < q max x < 3.8 for n = 4 and n = 6 respectively. 
Moreover, the saddle point value of the impact parameter, b s = 6 c (^6 c /n) _1// ( n+1 ) is only slightly 
less than b c for these parameter choices, failing to satisfy the condition b s <C b c needed for saddle 
point dominance of the phase integral, Eq. (|12p . These considerations tell us that, even with the 
maximal q consistent with use of the eikonal, the strong coupling, saddle point conditions are, at 
best, only marginally satisfied even at the highest, super GZK energies we consider. 

The sensitivity of the calculation to M$ and values is much more important to the estimate 
of event rates than the sensitivity to the use of the saddle point approximation, rather than the full 
EA. Since collisions that satisfy the EA conditions tend to be highly elastic, the consequences for 
event rate estimates based on the model assumptions are not severe, since the fraction of deposited 
energy is small, and the corresponding events are mixed with highly inelastic events at neutrino 
energies an order of magnitude or more lower, where the flux is presumably much larger. Thus the 
factor of ten difference in cross section between the Mq = 1 TeV and 2 TeV, or between the infinite 
and finite tension at fixed LSG is, in a sense, part of uncertainty in a "background" correction to 
inelastic processes, such as black hole formation, where the neutrino energy is nearly all converted 
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to visible energy. Only if black hole formation is severely suppressed, does the enhanced signal 
from elastic processes become a sensitive tool to detect LSG. 

In summary, our conclusions are that, in agreement with other studies in collider settings, the 
LSG effects are sensitively dependent on the choice of scale parameters - the gravity scale and the 
brane tension. The eikonal neutrino-nucleon cross section always rises sharply above the standard 
model versions, but the energy at which the onset occurs and the degree to which the new physics 
dominates, is strictly dictated by the input scales. In the regime we consider, 10 6 - 10 12 GeV, the 
saddle point dominance conditions are not well satisfied, and the approximation, though lending 
itself to simple expressions and intuitive pictures, noticeably underestimates the cross section. 
Having said that, the uncertainty in and Ms values is nevertheless the dominant effect. 
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S. is supported in part by NSF Grant No. PHY-0544278. D.W.M. receives support from DOE 
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Appendix 

The integral of Eq. (|12p can be expressed as a Mellin convolution: 

F n {y)=i r^rh(x)h(z)> (25) 

where 



o x \X 



A 0*0 = ^( e """ 1 ) (26) 

f 2 (x) = J (x). (27) 



We can calculate the Mellin transforms of Eqs. (|26p . (|27j) : 



f*( s ) = I" /ifoWdz = ~-r( — "\ e^( n - 1+ § ) (28) 
Jo n \ n J 



„ S -i^_ 2S - lr (§) 



/.(")- / h(x)x s ^dx= r(l _ j) - (29) 



o 



Upon defining f*(s) = /*(s)/|(s), the integral of Eq. ([2"5j) turns out to be a Mellin- Barnes integral: 

- iFn {y) = ^-l f*(s)y~ s ds (30) 



2ni _ 

that can be expressed explicitly in terms of Meijer G-functions. 
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The real part of Eq. (|30|) is 
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2vri 



1 f i2^r(f)r(^)sm[^(n-l + |)] 



2?ri 



71 



r(i-f) 



y s ds 



-y s ds, 



(31) 



2vr, 7nr(i-f)r(2=i + ^)r(i-^-^) 

where in the last step we have used Euler's duplication formula (sin(7rx) = 7r/[r(x)r(l — x)]). 

In the standard definition of Meijer G- functions [28] the coefficient multiplying s inside the 
Gamma functions has to be ±1. We therefore perform the substitution s — * 2ns to get 



vrr(n,)r (2 



2m J r(l - ns)T (2=1 + s) T (l 
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I) 



2« 



ds, 



(32) 



and use some known identities of the Gamma function to recast Eq. f|32|) into 
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ds. (33) 



The integral is now in the form of a Meijer G-function: 



Meijer G 
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A similar calculation for the imaginary part yields 
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In Eqs. (|34p and (|35p . the first line refers to the Mathematica input necessary to perform the 
calculation, while the second line refers to the standard mathematical notation of Eqs. (|14p and 
(El) M- 



13 



References 

[1] K. Greisen, Phys. Rev. Lett. 16, 748 (1966). 

G. T. Zatsepin and V. A. Kuzmin, JETP Lett. 4, 78 (1966) [Pisma Zh. Eksp. Teor. Fiz. 4, 114 
(1966)]; V. S. Berezinsky and G. T. Zatsepin, Phys. Lett. 28B (1969); Sov. J. Nucl. Phys. 
11, 111 (1970). 

[2] P. Biermann and G. Sigl, Lect. Notes Phys. 576, 1 (2001) |arXiv:astro-ph/0202425] . 

[3] W. Rhode et al. [Frejus Collaboration], Astropart. Phys. 4, 217 (1996); V. Balkanov et al. 
[BAIKAL Collaboration], |arXiv:astro-ph/0105269"} J. Ahrens et al, Phys. Rev. Lett. 90, 
251101 (2003) [arXiv:astro-ph/0303218l ; M. Ambrosio et al. [MACRO Collaboration], As- 
tropart. Phys. 19, 1 (2003) [arXiv:astro-ph/020318l| . 

[4] R. M. Baltrusaitis et al, Phys. Rev. D 31, 2192 (1985); S. Yoshida et al, Prepared for 27th In- 
ternational Cosmic Ray Conference (ICRC 2001), Hamburg, Germany, 7-15 Aug 2001; S. Hun- 
dertmark [AMANDA Collaboration], Prepared for 28th International Cosmic Ray Conferences 
(ICRC 2003), Tsukuba, Japan, 31 Jul - 7 Aug 2003; I. Kravchenko et al, Prepared for 28th 
International Cosmic Ray Conferences (ICRC 2003), Tsukuba, Japan, 31 Jul - 7 Aug 2003. 

[5] I. Kravchenko et al. [RICE Collaboration], Phys. Rev. D 73, 082002 (2006); S.W. Barwick 
et al. [ANITA Collaboration, Phys. Rev. Lett. 96, 171101 (2006); N. Letinen et al. [FORTE 
Collaboration], Phys. Rev. D 69, 013008 (2004); P. Gorham et al, [GLUE Collaboration], 
Phys. Rev. Lett. 93, 041101 (2004); "Limits on the Diffuse Neutrino Flux from the HiRes 
Data", W. Deng et al [HiRes Collaboration], Proceedings of the 29th International Cosmic 
Ray Conference. 

[6] A. Achterberg et al. [IceCube Collaboration], Phys. Rev. D 76: 027101 (2007), 507 (2004); A. 
Achterberg et al. [IceCube Collaboration], arXiv:astro-ph/0705.1315. 

[7] X. Bertou et al, Astropart. Phys. 17, 183 (2002). 

[8] H. Landsman [AURA Collaboration], Nucl. Phys. B Proceedings Suppl. 168, 268 (2007) 

[9] PGorham et al, Phys. Rev. D 72, 23002 (2005). 

[10] S. Bottai [EUSO Collaboration], Proceedings of the 27th International Cosmic Ray Conference, 
Aug. 2001, Hamburg, p. 848 (2001) 

[11] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 
larXiv:hep-ph/98033T5] . 



14 



[12] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [arXiv:hep-ph/9905221| . 

[13] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 257 
(1998) [arXiv:hep-ph/9804398] . 

[14] G. M. Frichter, D. W. McKay and J. P. Ralston, Phys. Rev. Lett. 74, 1508 (1995) [Erratum- 
ibid. 77, 4107 (1996)] |arXiv:hep-ph/9409433| . 

[15] R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Astropart. Phys. 5, 81 (1996) 
|arXiv:hep-p h/9512364 1 ; R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Phys. Rev. D 
58, 093009 (1998) |arXiv:hep-ph/9807264| . 

[16] P. Jain, S. Kar, D. W. McKay, S. Panda and J. P. Ralston, Phys. Rev. D 66, 065018 
(2002) |arXiv:hep-ph/0205052| ; S. Hussain and D. W. McKay, Phys. Rev. D 69, 085004 
(2004) |arXiv:hep-ph/0310091| ; S. Hussain and D. W. McKay, Phys. Lett. B 634, 130 (2006) 
[arXiv:hep-ph/0510083l . 

[17] M. Bando, T. Kugo, T. Noguchi and K. Yoshioka, Phys. Rev. Lett. 83, 3601 (1999). 
[arXiv:hep-ph/9906554^ ] . 

[18] M. Sjodahl and G. Gustafson, Eur. Phys. J. C53, 109 (2008) [arXiv:hep-ph/0608080] ; 
L. Lonnblad and M. Sjodahl, JHEP 0610, 088 (2006) |arXiv:hep-ph/0608210| . 

[19] J.J. Sakurai, "Modern Quantum Mechanics", Menlo Park, CA: The Benjamin/ Cummings 
Publishing Company (1985) 

[20] H. Cheng and T. T. Wu, "Expanding Protons: Scattering at High-Energies" Cambridge, USA: 
MIT-PR. (1987) 

[21] D. Kabat and M. Ortiz, Nucl. Phys. B 388, 570 (1992) |arXiv:hep-th/9203082| . 

[22] D. Kabat, Comments Nucl. Part. Phys. 20, 325 (1992) [arXiv:hep-th/9204103| . 

[23] G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B 544, 3 (1999) 
[arXiv:hep-ph/981129T] . 

[24] T. Han, J. D. Lykken and R. J. Zhang, Phys. Rev. D 59, 105006 (1999) |arXiv:hep-ph/9811350| . 
[25] R. Emparan, Phys. Rev. D 64, 024025 (2001) [arXiv:hep-th/0104009| . 

[26] G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B 630, 293 (2002) 
[arXiv:hep-ph/011216l] . 



15 



[27] R. Emparan, M. Masip and R. Rattazzi, Phys. Rev. D 65, 064023 (2002) [arXiv: 
hep-ph/0109287] . 

[28] A. Erdelyi et al., "The Bateman Manuscript Project", McGraw-Hill: New York, Toronto, 
London (1954), Vol. 1. 

[29] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000) |arXiv:hep-ph/9903282| . 

[30] S. Hannestad and G. G. Raffelt, Phys. Rev. Lett. 88, 071301 (2002) [arXiv:hep-ph/0110067] . 

[31] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 95, 161602 (2005) 
[arXiv:hep-ex/0506063] . 



16 



